J/i/jD*D* vertex from QCD sum rules. 



- 1—1 

X 



M.E. Bracco, M. Chiapparini 
Instituto de Fisica, Universidade do Estado do Rio de Janeiro, 
Rua Sao Francisco Xavier 524, 20550-900 Rio de Janeiro, RJ, Brazil 

F.S. Navarra and M. Nielsen 
Instituto de Fisica, Universidade de Sao Paulo, CP. 66318, 05389-970 Sao Paulo, SP, Brazil. 

We calculated the strong form factor and coupling constant for the J/ipD* D* vertex in a QCD sum 
rule calculation. We performed a double Borel sum rule for the three point correlation function of 
vertex considering both J ftp and D* mesons off-shell. The form factors obtained are very different, 
but they give the same coupling constant. 



PACS numbers: 14.40.Lb,14.40.Nd,12.38.Lg,11.55.Hx 



o 

o . 

An important topic of RHIC physics is charmonium production pj . The first results have already been reported 0] , 
| but they did not reach yet the quality of the previous lower energies measurements performed at CERN-SPS, where 
O ■ an anomalous suppression was found in the most central collisions Q. 

| /" ( | Understanding charmonium production and destruction in these highly complicated collisions is a challenge for 
theorists, and in the last six years an intense effort was dedicated to understand charmonium interactions with light 
hadrons. This subject can be divided into high (y/s > lOGeV) and low (yfs < lOGeV) energy interactions. The 
latter have been treated mostly with the help of effective Lagrangians. In this approach one has a better control of the 
relevant symmetries, which dictate the dynamics jj, la, 111 13 u| • However, very soon it was realized that working with 
£> | effective lagrangians requires a very good knowledge of the form factors associated with vertices involving charmed 
7 ~ I ■ mesons as, for example, the D*Dtt vertex, to mention the most famous one. Choosing a softer or harder form factor 
may change the final cross section up to two orders of magnitude. This dramatic change outshines the detailed 
discussion concerning the role played by gauge invariance, anomalous parity couplings and other dynamical features 
I ' of the interaction. 

, Until four years ago nothing was known about such form factors. At that time we launched a program of calculating 
these important quantities in the framework of QCD sum rules Ql. Since then we have been continuously attacking 
this problem and computing different vertices [ToL ITll IT2I fl3L Il4 ITU fill . In doing so, we have also improved the 
strategies to determine the coupling constants. One of them is experimentaly measured but the others have to 
be estimated using a vector meson dominance analysis of some measured decays, or using SU(4) symmetry relations. 
In some cases our results gave support to these more empyrical estimates and in some other they did not. 

As it frequently happens, during the execution of the project the original motivation was extended to a wider scope 
of questions concerning charmed mesons. A particularly interesting side product of our works, |l2l IT^j and also the 
present one, is the conclusion regarding the "resolving power" of a compact (J/ip) or a extended probe (/?) hitting the 
D meson. It was reassuring to observe how the J/ip behaves more like a point-like parton penetrating the D meson, 
whereas the p behaves more like a large hadron being able to measure the size of the D meson. All this information 
is encoded in the form factors. In what follows we will have a glimpse on how the J /ip "sees" a D* meson. From 
spin symmetry of the heavy quark effective theory (HQET) one would expect a similar behaviour as observed in the 
J/ipDD vertex [T^|. 

The present calculation is part of this project. Although the method is the same as used before, because there are 
three vector particles involved, the number of Lorentz structures is much bigger and the calculation is considerably 
more involved. On the other hand, since the three particles are heavy, we feel more confident about neglecting some 
higher dimension contributions to the operator product expansion (OPE). We will compare our results with the ones 
from other models 0, 0] and we will also check if the use of HQET symmetry is appropriate. 

Following the QCDSR formalism described in our previous works [ToL ITll Il2l|. we write the three-point function 
associated with the J/ipD* D* vertex, which is given by 

rl^W) = / d'xd'y e^e-^-p)-*< |T{jf (y)jf (0)}|0) (1) 

for an off-shell D* meson, and: 

r^fW) = / d 4 xd*y eV-'e-M-ti-y (0|T{jf *t(0)}|0) , (2) 



for an off-shell J/ip meson. The general expression for the vertices and @ has fourteen independent Lorentz 
structures. We can write T vafl in terms of the invariant amplitudes associated with each one of these structures in 
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the following form: 



IVab,??') = Ti(p 2 ,p' ,q 2 )g^p a +T 2 {p 2 ,p' -,q 2 )g^ a Pu + T 3 (p 2 ,p' , q 2 )g va p^ + T 4 (p 2 ,p' ,q 2 )g^p' a 

+F5(p 2 ,p' 2 ,q 2 )9v, a pl + F 6 (p 2 ,p' 2 ,q 2 )g ua p^ + r 7 (p 2 ,p' 2 ,q 2 )p^p„p a +T s (p 2 ,p' 2 ,q 2 )p^p' ly p a 

+^9(p 2 ,p' 2 ,q 2 )p^p'uP a + ri (p 2 \/ i \q 2 )p^Pup' a + Ti 1 (p 2 ,p' 2 ,q 2 )p^plp a + r 12 (p 2 ,p' 2 ,q 2 )p' fl p„p' a 
+^i3{p 2 ,p' 2 ,q 2 )PiipWa + T ii(p 2 ,p' 2 ,q 2 )p'f,plp' a (3) 

The correlator function, Eqs. Q an d can be calculated in two diferent ways: using quark degrees of freedom 
-the theoretical or QCD side- or using hadronic degrees of freedom -the phenomenological side. In the QCD side 
the correlators is evaluated by using the Wilson operator product expansion (OPE). The OPE incorporates the 
effects of the QCD vacuum through an infinite serie of condensates of increasing dimension. On the other hand, the 
representation in terms of hadronic degrees of freedom is responsable for the introduction of the form factors, decay 
constants and masses. Both representations are matched using the quark-hadron global duality. 

In the QCD side, each meson interpolating field appearing in Eqs. Q and (J2J) can be written in terms of the quark 
field operators in the following form: (x) — c(x)j v q(x) and jf(x) — c(x)^ fl c(x), where q and c are the up/down 
and charm quark field respectively. Each one of these currents have the same quantum numbers as the associated 
meson. 

For each one of the invariant amplitudes appearing in Eq.Q, we can write a double dispersion relation over the 
virtualities p 2 and p' , holding Q 2 — —q 2 fixed: 



T l (p 2 1 p' 2 ,Q 2 ) = -- / ds / 



du ■ 



pi(s,u,t 



(s — p 2 )(u — p' ) 



i = l, 



,14 



(4) 



where pi(s,u,Q 2 ) equals the double discontinuity of the amplitude Ti(p 2 ,p' ,Q 2 ), calculated using the Cutkosky's 



rule, and where s m i n = Am 2 in the case of D* off-shell, and s T: 



m; 



in the case of J/ ip off-shell, with m c being the 



mass of the charm quark. The invariant amplitudes receive contributions from all terms in the OPE. The first one of 
those contributions comes from the perturbative term and it is represented in Fig. ^ 
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FIG. 1: Perturbative diagrams for D* off shell (left) and J/ip off shell (right). 

In principle, we can work with any structure appearing in Eq.Q. But we must choose those which have less 
ambiguities in the QCD sum rules approach, which means, less influence of the condensates of higher dimension, 
and a better stability as a function of the Borel mass. A possibility could be work with the structures that present 
some symmetry in the momenta of the mesons, for example, Pp,p v p a or p'aP'vP'a- I n these cases, however, the quark 
condensate is absent. Since the quark condensate is the condensate of lower dimension we will choose structures in 
which the quark condensate contribute. These structures are g^ip' + p) a , g av (p' — p)^ and g afl (p' — p) v . 

In this paper we use the structure g^uip' +p)a, which presented a better behavior. The corresponding perturbative 
spectral densities which enters in Eq. (J1J is 



p( D *\ S ,u,Q 2 ) = 



s - t 



M + N 



2(C + D) + -(3m 2 - u - s) - 2G 



(5) 



for D* off-shell, and 



MM 



(s,u,Q 2 ) = - 



t — u — s 



U + A N -2(C + D) + ^(u + s-2m 2 ) + 2G 



(6) 
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for J/ip off-shell. Here A = A(s, u, t) = s 2 + t 2 + u 2 - 2st - 2su - 2tu, s = p 2 , u = p' 2 , t = -Q 2 and C, D, G, M and 
N are functions of (s,t,u), given by the following expressions: 

C = !^(l-cos^^cos^-^ (7) 

—3 

D = -^|L(l-cos 2 0)cos0 (8) 

G = -ix\k\ (1 - cos 2 0) (9) 

M = tK^IH) (10) 

N = 27T-1^-COS0 (11) 



m 

fit 



s + u — t , . 

Po = -^T7^- (12) 



l^l 2 - £ (13) 

2y/S 



-2 



-2 



|&| = fe -m 2 (15) 
cQs - = _ u + ^g-_2^ (16) 
2\f?\\k\ 

where m — and 77 = 1 for D* off-shell and m = m c and rj = —1 for J/^> off-shell. 

The contribution of the quark condensate which survives after the double Borel transform is represented in Fig. [3 
for the J I tjj off-shell case, and is given by 

r (J /V >) = m c(qq) 

(p 2 — m 2 )(p' 2 — m 2 ) 

where (qq) is the light quark condensate. For the D* off-shell there is a similar contribution, this time proportional 
to the mass of the light quark and to the cc condensate, both small, therefore, we will not consider it here. We also 
expect that the perturbative contribution be the dominant one in the OPE, because we are dealing with heavy quarks. 
For this reason, we do not include the gluon and quark-gluon condesates in the present work. 




D* ill) D* 

FIG. 2: Contribution of the qq condensate for J/ip off-shell case. 

The resulting vertex functions in the QCD side for the structure g^ip' +p)a are written as 

1 



so ds r du ^ 

4m? J m? 



for D* off-shell and 



rwi>)( v „'\ = L f ,, , , „„ 



ds I du 



(s — p 2 )(u — p' 2 ) 



2 (s-p 2 ){u-p' ) 



(18) 



(19) 
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for J/t/j off-shell, where, as usual, we have already transferee! the continuum contribution from the hadronic side to 
the QCD side, through the introduction of the continuum thresholds Sq and Uq |20|. 

The phenomenological side of the vertex function is obtained by considering the contribution of the J/ ip and two 
D* mesons in Eqs. and ©■ The effective Lagrangian for the J/ipD*D* interaction is given by p| 



■J/^D'D* — ig-ipD-D* 



-rd»Dt)D*» + D^^d^DZ - (d^)Dt)] , (20) 
from where one can extract the matrix element associated with the J/ipD*D* vertex: 
(2r(pi,a)D*(pa,i/)|J/V>(p3,/i)) = 

- -igj/^D'D* i(P2 -V\)^9av - {P3+P2)ng^ + <J>3 + Pi) n9 tia]^, /^iPs)^* D' (?!>£>* feO- ( 21 ) 

The meson decay constants, fo* and fj/^p, are defined by the matrix elements 

{0\f D ,\D*(p))=rn D ,f D ,e^(p) (22) 

and 

(0\j^m(p)) = mj^fj/^ip) (23) 

where e^,* and are the polarization vectors of the mesons D* and J/ip respectively. Saturating Eqs. QJ and @ 
with the J /ip and two D* states and using Eqs. (|2ip. (|22|l and l|23|l we arrive at 

v {j/i>) _ J-VVO (r >2\ fh*f.J/4> m D* m J/4: ( _ _ P»P»> 

K a » - .9,7/^*W^ ( p 2 + m 2 )0(Q2+TO2 ^ )(p ,2 + m ^ ) ( v 9* ^ 



_ P'uPv' \ _ _ gaga' 



(p + j/) a y V +{q-pfg°"^ ~ (p' ~ qT' 9 a ' u ' , (24) 



when the J/tp is off-shell, with a similar expression for the D* off-shell. The contractions of n' , v' and a' in the above 
equation will lead to the fourteen Lorentz structures appearing in Eq. (PJ). We can also see from Eq. (|24|1 that the 

form factor 9j/J^)» D , (Q 2 ) is the same for all the structures and thus can be extracted from sum rules written for any 
of these structures. 

The resulting phenomenological invariant amplitudes associated with the structure g^vip' +p) a are 



T {D * ] (n 2 r/ 2 <1 2 ) - n^ D "> (D 2 ) /p* /j/ifi M D' w J/^ 

ip,, (P ,P M ) 9j/^ D , D ,W )(pa + m i \. )(Q* + m l D .)(P* + m 2 D .) 



(25) 



for the D* off-shell, and 



y vh [P ,P ,W )-9 J/i)D * D *W ^pa +m ^j(Q 2+m 5 \ ){P' 2 +m%,) [ 1 



for J/ip off-shell. 

To improve the matching between the two sides of the sum rules we perform a double Borel transformation pof in 
the variables P 2 = — p 2 — > M 2 and P' 2 = —p' — ► M' 2 , on both invariant amplitudes T and r p h- Equating the results 
we get the final expression for the sum rule which allow us to obtain the form factors g^/]p D „ D *(Q 2 ) appearing in 
Eqs. H25fl - 126|l . where T is D* or J/ijj. In this work we use the following relations between the Borel masses M 2 and 
M f2 M^ = <L± for an off-shell D* and M 2 = M' 2 for an off-shell JM. 

The values of the parameters used in the present calculation are presented in TableQ] The continuum thresholds are 
given by sq = (m + A s ) 2 and uq = (m + A u ) 2 , where m is the mass of the incoming meson. We used the experimental 
value for the J/ip decay constant /m, and took fn* from ref. |2ll |. 

Using A s = A u = 0.6 GeV for the continuum thresholds and fixing Q 2 — 1 GeV 2 , we found a good stability of the 
sum rule for 9j/^d*d* ^ or ^ 2 * n ^ e interval 6 < M 2 < 18 GeV 2 , as can be seen in Fig. |3| 
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m c (GeV) 


m D .(GeV) 


mj/^i GeV) 


/WGeV) 


f J/ A GeV) 


{qq)(GeVf 


1.3 


2.01 


3.1 


0.240 


0.405 


(-0.23) 3 


TABLE I 


Parameters used in the calculation. 
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FIG. 3: Stability of g^ D , D , (Q* = lGeV), as a function of the Borel mass 



M 



In the case of g^//J^,* D * the interval for stability is 5 < M 2 < 10 GeV 2 , as can be seen in Fig. 0] 
Fixing A s = A„ = 0.6 GeV we calculated the momentum dependence of the form factors. We present the results in 
Fig. 03 where the circles corresponds to the gf/^ D , D *(Q 2 ) form factor (evaluated using M 2 — 9. GeV 2 ) in the interval 

where the sum rule is valid. The triangles are the result of the sum rule for the 9^'jj^» D * {Q 2 ) form factor (evaluated 
using M 2 = 6. GeV 2 ). 



In the case of the off-shell D* meson, our numerical results can be fitted by a monopolar parametrization, dashed 
line in Fig. in the following way: 



(m , 2) _ 6.04 



(27) 



As in Ref . [T^ . we define the coupling constant as the value of the form factor at Q 2 = —m^, where mj- is the mass 
of the off-shell meson. Therefore, using Q 2 = —m 2 D , in Eq (|27fl. the resulting coupling constant is: 



9,j/^d*d* = 6.75 



(28) 



In the case of the off-shell J/ip meson, our sum rule result can also be fitted by a monopole parametrization, which 
is represented by the solid line in Fig. [5j 



{Or) 



399.7 



Q 2 + 78.52 



Using Q 2 = —m 2 ^ in Eq (|29"|) we get: 



in a good agreement with the result in Eq. 



gj/ipD'D' — 5.80, 



(29) 



(30) 
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FIG. 4: Stability of g { j/^« D , (Q 2 = 1 GeV 2 ), as a function of the Borel mass. Here are represented the perturbative, condensate 
and total contributions. 
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FIG. 5: 9ji^£)*d* (circles and squares) and g\j/^, D » (triangles) form factors as a function of Q 2 from the QCDRS calculation 
of this work. The solid and dashed lines correspond to the monopole parametrization of the QCDSR data, using the two 
structures cited in the text, and the dotted line corresponds to the exponential parametrization. 



As said in the introduction, in principle we could work with any one of the structures appearing in Eq .@. To 
study the dependence of our results with the chosen structure, we show in Fig. 5, through the squares, the form 
factor gJJ^ D * D *{Q 2 ) obtained from the p^PvPa structure, using the same procedure described above. In this case the 
QCDSR results can be fited by an exponential parametrization, dotted line in Fig. 5, given by 



9 ( j%W(Q 2 ) = ^ Q/3 - 55 



(31) 
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which gives, at the D* pole, gj/^D-D* = 6.11. Although the parametrizations in Eqs. (|27|l and (|31|l are rather 
different, the behaviour of the curves is not so different, as can be seen in Fig. 5, and lead to consistent values for the 
coupling constant. These differences can be attributed to the uncertainties of the method. 

In order to study the dependence of our results with the continuum threshold, we vary A Si „ between 0.5 GeV < 
< 0.7 GeV in the parametrization corresponding to the J/ip off-shell case. As can be seen in Fig- EI this procedure 
gives us a uncertanty interval of 4.8 < gj/-^D*D* < 6.7 for the coupling constant. 




FIG. 6: Dependence of the form factor with the continuum threshold for the J/ip off-shell case. The lower curve (dash dot dot) 
corresponds to A s , u = 0.5 GeV and the upper one (dot) to A s ,u = 0.7 GeV. The dash correspond to A Sjli = 0.7, 0.6 GeV, the 
solid to A s ,„ = 0.6, 0.6 GeV, and the dash dot to A s ,„ = 0.6, 0.5 GeV 



Concluding, the two cases considered here, off-shell D* or J ftp, give compatible results for the coupling constant, 
evaluated using the QCDSR approach. Considering the uncertainties in the continuum thresholds we obtain: 

9 =6.2 ±0.9. (32) 

In Table ITTl we compare our result with those of other calculations. Our value is smaller than the values obtained 
using the constituent quark-meson model |l8j and a vector meson dominance approach [{| and is greater than the 
result obtained using a chiral model 0|. It is also in agreement, within the errors, with the QCDSR calculation for 
the J/ipDD coupling of Ref. 0, which is expected as a consequence of the spin symmetry of the HQET. 



Coupling 


This work 


Ref. [5] 


Ref. [14] 


Ref.flS] 


Ref. [191 


gj/i,D*D* 


6.2 ±0.9 


7.6 


5.8 ±0.8 


8.0 ±0.5 


4.9 



TABLE II: Values of the coupling using different approaches. Ref.[f| is a vector meson dominace calculation, Ref.|l4| is a 
QCDSR calculation in the same spirit as the calculation presented in this paper, Ref. |18|l uses the constituent quark model and 
Ref. uses a quiral model for computing the coupling. 

From the parametrizations in Eqs. i|28|) and l|29|l we can also extract the cutoff parameter, A, associated with the 
form factors. We get A ~ 2.2 GeV for an off-shell D* meson and A ~ 8.9 GeV for an off-shell J/ip. The cutoffs 
obtained here follow the same trend as observed in Refs. [ill IT2I IT3L IbH ] : the value of the cutoff is directly associated 
with the mass of the off-shell meson probing the vertex. The form factor is harder if the off-shell meson is heavier. 
Therefore, in the J/ipD*D* vertex, as in the J/ipDD, the J/ip behaves like a point- like parton penetrating the D* 
meson. 
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